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ABSTRACT 



, We have determined the physical mechanism responsible for the plasma in- 

stabilities, which was first found by Ramani and Laval (1978), associated with 
anisotropic velocity distributions induced by the temperature gradient in which 
there are growing low frequency transverse magnetic waves, even in the absence 
of background magnetic fields. We have shown that the physical mechanism re- 
sponsible for the growth of one of the modes is identical to the Weibel instability. 
The nonlinear saturation level of the instability is also provided by considering 
the wave-particle interactions. The non-linear evolutions of the magnetic fields 
after the saturation are speculated. The results are applied to the cold fronts 
which is one of the newly discovered structures in clusters of galaxies by the 
Chandra X-ray observatory. We predict the existence of the magnetic field of 
~ 10/uG tangential to the surface over the entire region of the cold front surface 
and that the heat conduction is significantly suppressed by the trapping of the 
electrons by the generated magnetic fields. The instability may provide a new 
possibility on the origin of cosmic magnetic field. 

Subject headings: galaxies: clusters: general — magnetic fields — conduction — 
inst abilit ies — plasmas 



1. INTRODUCTION 

The Chandra X-ray observatory revealed the bow-shaped discontinuities in the X-ray 
emitting hot plasmas in the clusters of galaxies (Markevitch et al. 2000). The plasma tem- 
perature sharply increases across the discontinuities, while the plasma density sharply de- 
creases across the discontinuities toward the same direction. These new structures are called 
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cold fronts. The cold fronts might be the contacting surfaces of two different origin plasmas. 
If the classical Spitzer heat conductivity is applied, the life time of the cold fronts should be 
~ 10 6 yr. Since it is much less than the expected ages of the cold fronts, that is > 10 8 yr, the 
heat conductivity in the cold fronts region must be significantly reduced from the classical 
Spitzer value (Ettori & Fabian 2000; Markevitch et al. 2000). Since there are thermal pres- 
sure jumps across the fronts (Vikhlinin, Markevitch & Murray 2001a), a transonic motion 
of the cold fronts is expected to maintain the pressure balance between the cold and hot re- 
gions. Vikhlinin et al. (2001b, 2002) pointed out that the magnetic field of 10/iG along the 
fronts should exist to explain the smoothness of the fronts otherwise the irregular structures 
are expected due to the Kelvin-Helmholtz (KH) instability. Since the required field strength 
is more than order of magnitude stronger than the intracluster magnetic field strength while 
the plasma density jump is at most factor of two across the cold fronts, the compression of 
the intracluster magnetic field as conserving the flux can not explain the existence of such a 
strong magnetic field in the cold fronts. 

A precise mathematical treatment of the thin hot plasmas based on the plasma kinetic 
theory found that the anisotropic electron velocity distribution induced by the temperature 
gradient drives the plasma instability, and it is postulated that the electron scattering by 
the excited plasma waves may significantly reduce the heat conductivity (Ramani & Laval 
1978). The low frequency transverse growing magnetic waves are excited even in the absence 
of an background magnetic field. Since the instability was found by Ramani and Laval 
(1978), the instability is referred to the Ramani-Laval (RL) instability in this paper. The 
applications of the RL instability to the cluster hot plasmas were examined by some authors; 
in cases of the plasma with background magnetic fields (Levinson & Eichler 1992) and the 
electron-ion two components plasma (Hattori & Umetsu 2000). Although the RL instability 
was said to be the Weibel like instability (Weibel 1959) since the growth of the waves is 
seeded by the anisotropic velocity distribution, there are plenty of the qualitative differences 
among the both instabilities as explained in section 3. Therefore, the both instabilities were 
distinguished (Gallev & Natanzon 1991) and it is very important to determine the physical 
mechanism of the RL instability. 

Since this paper is the first paper of our series of papers, we focus on the identification 
of the physical mechanism of the RL instability (in §2. and 4.) including the review of the 
Weibel instability (in Appendix) and the RL instability (in §3), the determination of the non- 
linear saturation level of the excited waves and the non-linear evolution of the excited modes 
(in §5), and the application to the cold fronts (in §6.). Section 7 is devoted to discussion. 



- 3- 



2. THE VELOCITY DISTRIBUTION FUNCTION IN THE PLASMA 

WITH THE HEAT FLUX 

In this section, how the anisotropic electron velocity distribution function is set up 
when the temperature inhomogeneity exists in the plasma, is discussed from physical point 
of view. The absence of a background magnetic field is assumed. Consider the temperature 
inhomogeneity in hot electron plasma with the temperature variation scale of L. Since the 
heat conduction carries the heat flux from the hotter to cooler regions, the heat flux along 
the temperature gradient takes a negative finite value; q oc (v\\v 2 f) < where the bracket 
denotes the integral over the velocity space and the subscript || denotes the component 
parallel to the temperature gradient. Therefore, the electron velocity distribution function 
must deviate from the Maxwell-Boltzmann f m = n (a;||)(7n; t h(£||) 2 )~ 3 / 2 exp(— v 2 /f t h(^||) 2 ); 
Af = f — f m 7^ 0. Here, no(x\\) is the electron number density, v t h(x\\) = ( ! 2kBT(x\\)/m e ) 1 ^ 2 
is the thermal velocity with the temperature T(x\\) and the electron mass m e . The above 
condition together with (Af) = 0, that is the number density conservation, and (v 2 Af) = 0, 
that is the energy conservation, restrict the form of the deviated part, Af, to be odd function 
of the velocity component along the temperature gradient, v\\. Of course, Af — > as 
\vu | — > oo. In Figure 1, the possible cases as Af are shown. From zero current condition, 
such as (v\\Af) = 0, the form like the type A is rejected since Af where v» > and v» < are 
respectively positive and negative, and then v\\Af > for all velocity and (v\\Af) > 0. On 
the other hand, the type B and C satisfy the zero current condition since Af curve acrosses 
the i>||-axis in each positive and negative v\\ region. The (v\\v 2 f) for the type B and C have 
finite values because this is the weighted mean of (v\\Af) weighted by v 2 which is larger 
when |f|j| is larger. The heat flux condition of (v\\v 2 f) < says that Af where v\\ — > — oo 
should be positive. Thus, the type C is only the possible form as the deviated part in the 
plasma where the finite heat flow exists. The relative amplitude of the deviated part to 
the Maxwell-Boltzmann should be e5x where e = A m f p /L and St = ST/T is the fractional 
temperature fluctuation, since the deviation may be induced by the Coulomb collision and 
the temperature fluctuation. 

The deviated part Af can be also deduced analytically. The Boltzmann equation is 

df , df e df 
at ox\\ m " ov\\ 

where E\\ is the zeroth electric field along the temperature gradient and the rhs is the Krook 
operator as the collision term and v = {k B T jm^) x l 2 j\ e is the Coulomb collision frequency 
with the Coulomb mean free path A e (Sarazin 1988). Hereafter, we describe the collision 
frequency as v = fth/A m f P with the parameters i> t h = (2/csT/m e ) 1//2 and A m f p = V2\ e - For 
simplicity, the pressure balance is assumed. Then, E\\ = (Ramani & Laval 1978). If the 
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perturbative treatment is applicable to describe the system, the distribution function could 
be expanded in e8t (Chapman & Cowling 1960) as, 



where (j = 1,2---) describe the deviation of the distribution function from the Maxwell- 
Boltzmann in order of (e&r) 3 '. This expansion is known as the Chapman- Enskog expansion. 
Therefore, the electron distribution function up to the first order in e5r is obtained as 



The form of the deviated part is essentially the same as the type C shown in Figure 1. We 
would like to note that adopting the Krook operator is not essential for determining the form 
of the deviated part. 



The diagnostics of the Ramani-Laval (RL) instability are summarized contrasting with 
the Weibel instability. As summarized in the Appendix, the temperature anisotropy excites 
the transverse magnetic waves due to the Weibel instability. The excited wave is standing 
wave with zero phase velocity. The amplitude of the wave grows when the temperature 
perpendicular to the wave vector is higher than the temperature parallel to the wave vector. 
The amplitude of the wave is damped in the opposite situation. 

Ramani and Laval (1978) studied the stability of the plasmas when the temperature 
distribution is not homogeneous. The non-equilibrium electron velocity distribution function 
is deduced from the assumed temperature distribution self consistently using the Chapman- 
Enskog expansion. This is one of the prominent difference from the Weibel case where the 
anisotropic electron velocity distribution function due to the temperature anisotropy is given 
by hand. The difference of the deduced non-equilibrium velocity distribution function from 
the equilibrium one, that is the Maxwell-Boltzmann distribution, is skewed in the direction 
of the temperature gradient and is the odd function of the velocity component along the 
temperature gradient. The velocity distribution functions in the both cases are anisotropic. 
However, the skewed nature found in the RL case is not found in the Weibel. Although the 
anisotropic velocity dispersion is essential for the Weibel instability, the velocity dispersion is 
isotropic in the RL case. Ramani and Laval (1978) performed the linear stability analysis of 
the plasmas with the deduced non-equilibrium velocity distribution following the procedure of 
the plasma kinetic theory. Two independent modes appear. The two modes are distinguished 
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whether the magnetic field is in or not in the plane made by the temperature gradient and 
the wave vector, where the first one is named mode 2 and the second one is named mode 1. 
The dispersion relation of the both modes have a real part as 

Ld r = ^-/o> t hCOS#. 

The imaginary part of the mode 1 is obtained as, 

e 2 5 2 1 / c \ 2 

7 = — ^kv th (3cos 2 9 - 2sin 2 #) - — — k 3 v th , 
80r 0r \uj p J 



and the imaginary part of the mode 2 is obtained as, 

3 ^T, 2 Q 1 ( C V,3 

7 = -s-^KUthCos 9 = — k vth, 



where u; p is the electron plasma frequency, c is the speed of light and 9 is the angle between 
the direction of the temperature gradient and the wave vector. The existence of the real 
part in the dispersion relation of the RL instability is one of the difference from the Weibel 
case. Since the phase velocity along the temperature gradient uj r /kcos9 = ^v th only takes 
a positive value, the wave propagates only one way from the low temperature region to the 
high temperature region. Therefore, the excited waves do not carry the heat from the hot to 
the cold region. Although the mode 2 is a pure transverse wave, the longitudinal component 
of the electric fields, E k , has a non-zero value for the mode 1 as 

E k ~ —sm9 — B z . 

4 c 

This is an another prominent difference from the Weibel case. The wave grows the most 
rapidly when the wave vector is parallel to the temperature gradient. The wave is damped 
when the direction of the wave propagation is perpendicular to the temperature gradient. 
The dependence on the k of the growth rate is the same as that of the Weibel. 

In a sence that the anisotropic velocity distribution is the driving force of the insta- 
bility, the both instabilities are similar and therefore the RL instability was said to be 
the Weibel like. However, there are plenty of the qualitative differences among the both 
instabilities as explained in above. Therefore, the both instabilities were distinguished 
(Gallev & Natanzon 1991) and we try to identify the physical mechanism of the RL in- 
stability in the next section. 



4. THE PHYSICAL MECHANISM OF THE RL INSTABILITY 



The physical mechanism of the RL instability can be understood from the nature of the 
velocity distribution function when the finite heat flow exists, as illustrated in Figure 2. In 
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the following discussion, strictness of the numerical factor is ignored. The peak position is 



velocity of any low frequency magnetic transverse wave must be close to the velocity at where 
the distribution function has a peak value otherwise a finite net electric current is induced by 
the wave magnetic fields. This explains why the excited waves by the RL instability have the 
phase velocity of eSrVth, and travel one way direction rising up the temperature gradient and 
the waves can not travel in the direction of the perpendicular to the temperature gradient. 
The reason why the excited waves by the Weibel instability are standing wave, is simply 
because there is no shift in peak position of the velocity distribution function in the Weibel 
case. The peak value is increased by ~ 1 + (e5 T ) 2 from that of the pure Maxwell-Boltzmann 
case since the amplitude of the deviated part takes a value of e5rv\\ ~ (e<5r) 2 fth relative to the 
Maxwell-Boltzmann and the value of the Maxwell-Boltzmann part is nearly same as the peak 
value at v\\ ~ e^-r^th *C fth- Since the total electron number density must be unchanged, for 
the observer comoving with the waves this can be interpreted as the decrease of the effective 
electron temperature in the direction of the temperature gradient, TV by ~ 1 — (e<5r) 2 relative 
to the temperature perpendicular to the temperature gradient, T± (Fig. 2). The growth of the 
magnetic waves in the RL instability is, therefore, due to essentially the same mechanism 
as the Weibel instability (Weibel 1959; Fried 1959) in which the temperature anisotropy 
is the driving force of the instability. Consider the waves traveling nearly parallel to the 
temperature gradient. In this case, T^ ^ ~ Tj_ and T,, g ~ T| where T ± g and Tj, £ are the 
temperature components perpendicular and parallel to the wave vector for the observers 
comoving with the waves, respectively. Since T ± g > Tj, g, the waves can grow. As a result, 
the direction of the magnetic field generated by the instability is almost perpendicular to the 
temperature gradient. The growth rate of the mode which travels toward the direction of 
the temperature gradient with wavenumber k is obtained from that of the Weibel instability 



The growth rate gets the maximum value of 7 max ~ (e<5T) 3 (fth/cV P at k — /c max ~ €5tuj p /c. 
When the direction of the wave vector is perpendicular to the temperature gradient, T ± g = 
Tj| and Tj, g = T±. Since T ± £ < Tj, £ in this case, the wave can not grow. These results are 
exactly the same as the results led by the plasma kinetic theory except numerical factors 
(Ramani & Laval 1978; Hattori & Umetsu 2000; Okabe & Hattori in preparation). 



shifted toward positive v\\ direction and the amount of the shift is v 



tSrVth- The phase 



(Krall 1973), 
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The above discussion gives a complete physical explanation for the mode 2 of the RL in- 
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stability. However, these are not satisfactory to the mode 1 in which the longitudinal electric 
field component appears. Unfortunately, we have not yet gotten the physical explanation 
for this. This is still an open question. 

5. ON THE NONLINEAR EVOLUTIONS OF THE RL INSTABILITY 

The nonlinear saturation level of the excited wave is estimated assuming that the wave- 
particle interaction determines the saturation level. The fundamentals are illustrated in 
Figure 3 (Ramani & Laval 1978; Gallev & Natanzon 1991). Once the Larmor radius of an 
electron gets shorter than the wavelength of the growing mode, the electron is trapped by the 
magnetic field of the wave and the magnetic flux enclosed by its orbit becomes finite. Then, 
the kinetic energy of the trapped electron starts monotonically increase as the growth of the 
magnetic field strength, since the increase of the magnetic flux enclosed by the electron orbit 
cause the induction electric fields which accelerate the electron as the betatron accelerator. 
Once the Larmor radius of typical thermal electrons, r L ~ v^u^ 1 , gets shorter than the 
wavelength of the fastest growing mode, that is ri/c max < 1, the increase of the kinetic 
energy of the electron system becomes significant if the waves still continue to grow. Since 
this finally violates the energy conservation, the growth of the magnetic field strength must 
be saturated when r L k max ~ 1. 

The evolution of the magnetic fields after the nonlinear saturation could be described 
as follows. Some numerical simulations which follow the evolution of the Weibel instability, 
showed that the strength of the magnetic field driven by the Weibel instability decreases 
after it gets the maximum value (Morse & Nielson 1971). This can be understood as fol- 
lows. As the magnetic field grows, the electron velocity distribution is isotropitized and it 
becomes difficult to maintain the electric current field which supports the magnetic field 
of the waves. In these simulations, the system is assumed to be isolated and the initial 
anisotropic velocity distribution function let be free to evolve to isotropic one. On the other 
hand, in the case of the RL instability there is a driving force which maintains the anisotropy 
of the velocity distribution function. As far as the temperature gradient is not disappeared, 
the finite heat flux transports the heat from the hot to the cold region and the anisotropic 
velocity distribution function discussed in section 2 is maintained. Therefore, the decrease 
of the magnetic field strength after the nonlinear saturation as found in the Weibel instabil- 
ity may not be occurred in the RL instability. It is expected that the generated magnetic 
field strength is kept to be the saturated value for the life time of the temperature gradient. 
Wallace and Epperlein (1991) performed the numerical simulations to follow the evolution 
of the Weibel instability when an initial anisotropic distribution function is maintained by 
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an external source. They showed that the magnetic field strength is kept constant value 
after the saturation. Their results support the above expectation for the RL instability. 
There are several indicative numerical simulations concerning the organization of the glob- 
ally connected magnetic fields from the wavy fields generated by the instability. The wavy 
magnetic fields generated by the Weibel instability evolve into the longer wavelength modes 
after the saturation (Lee 1973; Sentoku et al. 2000; Sentoku et al. 2002). This result indi- 
cates that the excited wavy magnetic field automatically evolves into the globally connected 
fields. This could be understood as follows. After the magnetic field strength gets the satu- 
rated level, the electric current field starts to play as an individual electric beams every half 
wave length (Fig. 4). Each beam is surrounded by the azimuthal magnetic field generated 
by the current beam itself. The electric beams interact each other via the Ampere's force 
(Sentoku et al. 2000; Sentoku et al. 2002). The beams directed to the same direction are 
attracted each other and automatically gather. Finally, they merge into larger one beam. 
Since the physical mechanism of the growth of the RL instability is the same as the Weibel 
instability as shown in section 4, the same evolution is expected even in the RL case. 

Although the reduction of the heat conductivity was originally considered due to the 
electrons scattering by the waves generated by the RL instability (Ramani & Laval 1978; 
Hattori & Umetsu 2000), this may not be the case. As discussed in above, the wavy magnetic 
field generated by the instability could tend to form the global magnetic field automatically. 
Therefore, the suppression of the heat conductivity may be determined by the trapping 
of the electrons by the organized magnetic field. To estimate the suppression of the heat 
conductivity quantitatively, we have to know what is the final structure of the magnetic 
field due to the self-organization. The detail nonlinear studies, numerical simulations for 
example, are desired to answer the questions. 



6. APPLICATION TO THE COLD FRONT IN A3667 

The above results are applied to the cold front found in the cluster of galaxies A3667. 
The temperature and density of the electron changes from ksT c = 4.1 ± 0.2 keV and n ejC = 
3.2 ± 0.5 x 10~ 3 cm" 3 to k B T h = 7.7 ± 0.8 keV and n e , h = 0.82 ± 0.12 x 10~ 3 cm' 3 across 
the cold front (Vikhlinin, Markevitch & Murray 2001a). The width of cold front is L ~ 
5 kpc at maximum. The Coulomb mean free path of the thermal electron in the cold 
front region is A e = 5.4(T/T avc ) 2 (n e /n avc ) _1 kpc where T ave = (T h + T c )/2 and n ave = 
(^ e ,/i + n e>c )/2. Thus, e5x ~ 1 in the cold front of A3367. The growth time scale of the 
unstable mode in the cold front is 7~ ax ~ O.lsec. Although the exact time scale over which 
the instability can generate magnetic fields with the saturation values over the scale of the 
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interfaces is difficult to estimate, that must be several tens times longer than the growth time 
scale. Therefore, it is expected that the magnetic field is generated almost instantaneously 
compared with any other dynamical time scales in the cluster of galaxies. By applying the 
saturation level discussed in the previous section, the saturated magnetic field strength in 
the cold front should be B SP; ± ~ ^JnyJuekBTedr — 8 (T/T avc ) 3 / 2 (n/n avc ) -1 / 2 /iG where the 
exact dispersion relation obtained by Ramani and Laval (1978) was used. The obtained 
value agrees surprisingly well with the speculated value of 10 /iG based on the stability 
consideration against the KH instability. The direction of the generated magnetic field 
should be almost parallel to the cold front surface due to the characteristics of the mode. 
As discussed in section 5, the generated wavy magnetic fields could tend to evolve into the 
globally connected field. Therefore, we expect that the cold front surface is covered by the 
globally connected magnetic fields directed tangential to the surface with the strength of 
10 /iG and the KH instability is suppressed by these fields. Although to suppress the KH 
instability the existence of the magnetic fields are required only in the tail of the cold fronts 
(Vikhlinin, Markevitch & Murray 2001b), our model predicts the existence of such a strong 
fields all over the cold front where the temperature jump exists. The generated magnetic 
field may significantly reduce the electron mean free path in the direction of the temperature 
gradient and could suppress the heat conduction all over the cold fronts. 



7. DISCUSSION 

We have successfully determined the physical mechanism responsible for one of the two 
independent modes of the RL instability, and shown that the growth mechanism is identical 
to the Weibel instability which is well-known as one of the generation mechanism of the 
magnetic field. Therefore, the RL instability can be also considered as the generation mech- 
anism of the magnetic fields in the astronomical situations. The nonlinear saturation level 
of the instability is estimated by considering the wave-particle interaction. The evolutions 
of the magnetic fields after the saturation are speculated by referring the previous numerical 
simulations which followed the non-linear evolution of the Weibel instability. The generated 
fields might be self-organized and evolve into the globally connected magnetic field. The 
results are applied to the cold front in A3667. The existence of the magnetic field of ~ 10/iG 
tangential to the surfaces in all over the cold front is predicted. The results surprisingly agree 
with the predicted nature of the magnetic field in the cold front to explain the smoothness 
of the cold front surface structure avoiding the KH instability. Although the suppression of 
the KH instability requires the existence of the strong magnetic field only in the tail of the 
cold front, our model predicts the existence of such a strong fields all over the cold front 
where the temperature jump exists. The significant suppression of the heat conduction in 
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the cold front is also expected by the trapping of the electrons by the generated magnetic 
fields. Thus, our model predicts that the bow-shaped structure of the cold fronts should not 
partially disappear and the whole structure should be maintained. 

Of course further studies on the RL instability are desired. The non-linear studies using 
the numerical simulations specificated to the RL instability and the cold front are one of 
the most important remaining subjects. The identification of the physical mechanism of 
the mode 1 how the longitudinal electric field is generated, is another important remained 
subject. 

The physics of the RL instability discussed in this paper should be universal and the 
RL instability may play an important role in various hot thin plasmas in the universe. For 
example, the temperature inhomogeneities of e5x ~ 0.04 distributed in whole cluster has 
been found in almost all clusters of galaxies. The RL instability predicts the existence 
of the cluster halo magnetic field with strength of B ~ 0.04 x y/n e k B T ~ 0.3/iG. Sur- 
prisingly the predicted value just coincides with the observationally reported cluster halo 
magnetic field strength. The details of this topic are reported in the force coming paper 
(Okabe & Hattori in preparation). We believe that the RL instability may provide a new 
fruitful possibility on the origin of the cosmic magnetic field. 
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A. Appendix : THE WEIBEL INSTABILITY 

The Weibel instability is driven by an anisotropic temperature distribution and generates 
the growing transverse magnetic standing waves (Weibel 1959). This instability is well- 
known as the mechanism of the magnetic fields generation from zero initial magnetic field. 
Consider the plasma with an anisotropic temperature such as 



/ = n 



■ cxp 



J th, 



7 th,J 



where the subscripts ||,_L respectively denote the parallel and perpendicular direction to 
the wave vector. The dispersion relation is obtained by the linear plasma kinetic theory 
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(Krall 1973) as 



LU r = 0, 




The waves has no real part of the frequency. The waves grow when T± > Ty. The excited 
waves are transverse magnetic waves with no electric field. 

The physical mechanism of the Weibel instability is first identified by Fried (1959). 
Consider the simple situations as shown in Fig. 4. The electrons and their initial velocities 
are expressed as the filled circles and the dotted arrows. In the case A, the simple anisotropic 
initial electron distribution function 

fo(v) = n u5(v\\)5(v 2 ± - u 2 )5(v z ) 

is assumed where Uq is the electron number density. This represents an extreme case of 
T± 7^ 0,Tj| = 0. In the case B, the anisotropic initial electron distribution function 

f (v) = n w5(vf - w 2 )S(v±)S(v z ) 

is assumed. This represents an extreme case of T± = 0, Tji ^ 0. The directions of the 
perturbed magnetic fields are perpendicular to the paper : & <g> represent that the fields 
project from the paper and rush in the paper, respectively. The size of the circle is propor- 
tional to the field strength at each position. 

Firstly, the qualitative explanation of the physical mechanism is reviewed. The electrons' 
orbits are deflected due to the introduction of the perturbed magnetic field as shown by the 
dashed arrows. The amount of the deflection is proportional to the field strengths at the 
each position. We first discuss the Case A. Consider the electrons (a) and (a') which initially 
have the negative velocity in the x± direction. They respectively carry the positive electric 
currents in the x± direction out from and into the shaded region. Because of the difference of 
the field strengths, the net positive current in the x± direction is carried out from the region 
in a certain infinitesimal time interval. In other words the net negative electric current 
is carried into the region by the perturbed motion of these electrons. The electrons (b) 
and (b') also carry the net negative electric current into the shaded region. Therefore, the 
net negative electric current is developed in the shaded region due to the injection of the 
perturbed magnetic field. Since the values of the induced electric currents are larger where 
the gradient of the magnetic field strength is larger, the current fields illustrated in Fig. 4 
are set up where the length of the arrows are proportional to the amplitudes of the currents. 
The magnetic field is produced around each electric current according to the Ampere's law 



- 12 - 



as illustrated in Fig.4. The growth of the magnetic field is determined by the superposition 
of these fields. For example, at point X\ the net magnetic field generated by the current 
fields is normal to the paper toward readers and amplifies the injected perturbed field. Since 
the excited fields amplify the injected perturbed fields everywhere, the perturbed field grows 
in the Case A. On the other hand, the net induced electric currents direct opposite direction 
to Case A in the Case B. Therefore, the perturbed magnetic field is damped in the Case 
B. Hence, when u > w (T± > Ty) the magnetic field perturbation grows and when u < w 
(T_j_ < Tj|) the magnetic field perturbation is damped. While this simple explanation shows 
qualitatively how the growth occurs, the role of the electromagnetic induction for the growth 
of the instability, which was neglected in the above discussion, should be examined since it 
acts to reduce the growth of the magnetic field. 

Next, the semi quantitative explanation of the instability is given to see the role of the 
induction. For this purpose, only the case A is considered. Suppose that an initial perturbed 
magnetic field is 



B z = Be 



ikxn 



The x\\ and x±_ components of the equation of motion of the electron are given by 

dv\\ v± 
m e — = -e — B z , 
at c 

dv± v\\ ie dB z 

Tn P —r— = e — B z — 



dt c kc dt ' 

respectively, where the last term in x± component represents the force due to the inductive 
electric field. Therefore, the time derivative of the electric current flux is given by 

9 < J-L v \\ > _ „2 u 2 r 
at cm e 

where j± = —ev± is the x± component of the electric current carried by an single electron 
and < X >= — f d 3 vX(v)f(v) is the average over the velocity space. This current flux and 
the inductive electric field cause a change in the mean value of j± as 

d < j± > d < j ± v\\ > 2 • 1 dB z 



dt dx\\ m e kc dt 



The Ampere's law provides 



ick 

n <J±> = —r-B g , 

47T 
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where the displacement current is neglected since the low frequency mode is only considered. 
By combining these equations, we obtain 

Ji 

1 H — 

c 2 k 2 

The dispersion relation is obtained as 




It shows that the perturbation can grow. The second term in the [ ] of the left hand side 
is coming from the induction term. Therefore, the electromagnetic induction actually plays 
a role of the inertia for the growth of the perturbation and reduces the growth rate, but 
never stop the growth of the perturbation. The obtained growth rate is exactly the same as 
the result deduced by the full linear analysis based on the plasma kinetic theory where the 
induction effect is also taken into account (e.g. Melrose 1986). 



d 2 B z 
dt 2 



^u 2 B z 
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Fig. 1. — The possible candidates as the deviated part of the velocity distribution function 
from the Maxwell-Boltzmann. 
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Fig. 2. — The solid line is the v\\ section of the total velocity distribution function f m + in 
the plasma with the temperature gradient. The peak position is shifted by ~ e^r^th from the 
Maxwell-Boltzmann. The peak value is increased by ~ 1 + e 2 5^ compared with the Maxwell- 
Boltzmann. For comparison, the Maxwell-Boltzmann velocity distribution function shifted 
by ~ e^j-fth drawn in the dashed line. The velocity distribution function gets thinner in vn 
direction. This can be interpreted as the decrease of the effective temperature by ~ 1 — e 2 5^ 
in the direction of the temperature gradient. 
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when r L > k max 
B 




Fig. 3. — The nonlinear saturation by the wave-particle interaction. The top panel: When 
t l > k m l x , the thermal electrons travel throughout the waves but their orbits are randomly 
disturbed by the wavy magnetic fields. The bottom panel: Once tl becomes smaller than 
^max> the thermal electrons are trapped by the fields and feels net non zero fields. 
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Fig. 4. — The physical mechanism of the Weibel instability 



